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We construct series expansions for the scaling variables (which transform multiplicatively un- 
der a renormalization group (RG) transformation) in examples where the RG flows, going from an 
unstable (Wilson's) fixed point to a stable (high-temperature) fixed point, can be calculated nu- 
merically. The examples are Dyson's hierarchical model and a simplified version of it. We provide 
numerical evidence that the scaling variables about the two fixed points have overlapping domain of 
convergence. We show how quantities such as the magnetic susceptibility can be expressed in terms 
of these variables. This procedure provide accurate analytical expressions both in the critical and 
high-temperature region. 
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I. INTRODUCTION 

In many statistical mechanics or field theory problems one needs to calculate the macroscopic features of a system 

■ in terms of some microscopic parameters. In field theory language, the microscopic parameters are called the bare 
parameters and the macroscopic behavior is encoded in the low momentum n-point functions which can be used to 

T-H ■ define the renormalized quantities. Expressing the renormalized quantities in terms of the bare ones is a notoriously 
K*" i difficult problem. 

The renormalization group (RG) method §,| was designed in part to tackle this problem. However, its practical 
^0 [ implementation for realistic lattice models such as spin models or lattice gauge theories is still a formidable technical 

■ enterprise. Near RG fixed points, expansions are often available. The main problem consists in interpolating between 
' fixed points. An interesting attempt to model the interpolation between the weak and strong coupling regime of 

an asymptotically free gauge theory can be found in a recent publication of the QCD-TARO collaboration H. The 
interpolation was performed using the Monte Carlo simulation following a method developed in Ref . 0| . One delicate 
J step of this method is the choice of a finite number of operators which can be used to express the renormalized action 
^ ' after successive RG transformations. 

I ' In this article, we address the question of the interpolation between fixed points in models where the block spin 
Oh! method closes exactly and where the RG transformation can be calculated easily by numerical methods. One model 
f-H satisfying these requirements is Dyson's hierarchical model |^J^, for which only the local measure is renormalized 
• • I according to a rather simple integral equation. This integral equation, briefly reviewed in section |^ can be treated 
. ^ ' with numerical integration methods |7|^ , or by algebraic methods |^ based on the Fourier transform of the original 
K> . equation. The second method leads to very accurate determinations of the critical exponents JTo| , pT| . 
^ ' One could say that the hierarchical model (HM) is "numerically solvable", in the sense that for a given set of 
^ bare parameters, one can calculate the zero-momentum n-point functions. However, one has to repeat the tedious 
numerical procedure for every particular calculation and analytical expressions are more desirable. In this article, we 
show that this goal can be achieved by constructing the scaling variables corresponding to the fixed points governing 
the RG flows in the high-temperature (HT) phase. 

We use the terminology "scaling variables" for quantities which transform multiplicatively under a RG transfor- 
mation. In Ref. where these quantitites were introduced, Wegner call them "scaling fields". Indeed, in order to 
emphasize that they are functions of the bare parameters, we would prefer to call them "scaling functions" . However 
since this term is already reserved for other quantities, we will use the terminology "scaling variables" which seems 
to have prevailed over the years. 

The scaling variables can in principle be constructed in the vicinity of any fixed point. Near a given fixed point, 
they are simply the eigen-directions of the linearized RG transformation. When moving away from the fixed point, 
non-linear terms need to be added order by order to maintain a multiplicative renormalization. What is the domain 
of convergence of this expansion? Can we combine two expansions in order to follow the flows in a crossover region? 
These questions are in general very difficult to answer, but if the domains of convergence of the expansions fail 
to overlap, the whole approach is useless. The empirical results that we will present indicate that the domains of 
convergence of various expansions do overlap. These results should be considered as an encouragement to pursue this 
approach for other models. 
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To limit the technical difficulties of constructing these functions to orders large enough to get an idea about the 
asymptotic behavior of their series expansions, we have imposed some limitations on the calculations presented below. 
First, we only consider the HT phase. Second, we mostly consider the flows from the unstable fixed point to the 
HT fixed starting along the unstable direction. The perturbative relaxation of this second condition will be briefiy 
discussed. Third, we only discuss the magnetic susceptibility and not the higher moments. 

The article is organized as follows. In section ||, we review the basic facts about the HM and its RG transformation. 
We present a variant of the truncation method proposed in Ref. which will be used in the rest of the article. We 
also clarify the relationship between the truncation and the HT expansion. 

In the next four sections, we present, explain and illustrate the main ideas with a simplified one-variable model 
iQ where all the calculations are not too difficult. This model is simply a quadratic map with two fixed points, one 
stable and one unstable. The one- variable model is presented and motivated in section III, the scaling variables are 
constructed in section IV and their convergence studied in section We want to make clear that all our assertions 
concerning the convergence of series are based on the analysis of numerical values of the coefficients rather than on 
analytical results concerning these coefficients. The susceptibility is calculated in section Vl. The most important 
result for these four sections is illustrated Fig. |ll| which shows that scaling variables can be constructed accurately 
in overlapping regions. 

The rest of the article is devoted to generalizing the constructio n for the HM, for which the RG transformation 
can be approximated by a multivariable qu adrati c map. In section VH, we show how to choose the coordinates in 
order to solve the linear problem. In section VIII, we present approximations which allow one to calculate nonlinear 
expansions for the scaling variables, 
important result is illustrated in Fig. 



Finally, the questions of convergence are discussed in section IX 
P2| which indicates overlapping domains of convergence. 
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In more general terms, our article addresses the question of the non-linear behavior of the RG flows in models 
where non-linear expansions are calculable and can be compared with numerical solutions. In the first approximation, 
we have a linear theory. Trying to go beyond this first approximation has some similarities with trying to make 
perturbation near integrable systems in classical mechanics. 



II. DYSON'S HIERARCHICAL MODEL 



In this section, we review the basic facts about the RG transformation of the HM to be used in the rest of the 
paper. In order to avoid useless repetitions, we will refer the reader to Ref. § for a more complete discussion. In the 
following, we will emphasize new material such as the various possibilities available for the truncation procedure and 
the relationship between the truncation and the HT expansion. 



A. The RG transformation 



The energy density (or action in the field theory language) of the HM has two parts. One part is non-local (the 
"kinetic term") and invariant under a RG transformation. Its explicit form can be found, for instance, in section II of 
Ref. The other part is a sum of local potentials given in terms of a unique function V{(p). The exponential e"^*^"^^ 
will be called the local measure and denoted Wo (</>)■ For instance, for Landau-Ginsburg models, the measures are of 
the form Wo((/)) = e"^"^ "^"^ , but we can also consider limiting cases such as a Ising measure Wo{4>) = 5{4>'^ ~ !)• 
Under a block spin transformation which integrates the spin variables in "boxes" with two sites, keeping their sum 
constant, the local measure transforms according to the intergral formula 

W^^M = '^eWm^l^r-^^' J dc^'W^ (^) W„ (^) , (2.1) 

where (3 is the inverse temperature (or the coefficient in front of the kinetic term) and C„+i is a normalization factor 
to be fixed at our convenience. 
We use the Fourier transform 

/dk 
^e^'=*W„(fc) . (2.2) 

We introduce a rescaling of fc by a factor u/ ^ where u and s are constants to be fixed at our convenience, by defining 

i?„(fc) ^ wA , (2.3) 
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In the following, we will use s = Ij \fc. For c = 2^~^/^, this corresponds to the scaling of a massless gaussian field in 
D dimensions. Contrarily to what we have done in the past, we will here absorb the temperature in the measure by 
setting u — \f^. With these choices, the RG transformation reads 



= C'n+i exp 



1 9^ 



Rn 



V2k 



(2.4) 



We fix the normalization constant C„ so that i?,i(0) — 1. For an Ising measure, Ro{k) = cos{^/]3k), while in general, 
we have to numerically integrate to determine the coefficients of Ro{k) expanded in terms of k. 
If we Taylor expand about the origin. 



1=0 



(where a„^o — 1) then the finite-volume susceptibility is 



(2.5) 



- - --f (! 



The susceptibility x is defined as 



X'- 



lim 



Xn 



(2.6) 



(2.7) 



The susceptibility tends to a finite limit for < /3 < /3c, where (3c is a constant depending on c. For (3 equal to or 
larger than (3c, the definition of x requires a subtraction (see e. g., Ref. [ p5[ for a practical implementation). In the 
following, we will only consider the HT phase {(3 < (3c)- 

We have eliminated (3 from the mapping , moving it to the initial local measure. The mapping has then fixed points 
independent of the temperature. One of them is the "universal function" U{k) found empirically in Ref. and 
which can be written with great accuracy using numerical coefficients provided by Koch and Wittwer in Ref. [ p^ . 
The initial measure as a function of (3 determines a line of initial conditions in the parameter space, running from the 
HT fixed point, where all the initial parameters are zero, to the critical surface {(3 = (3c)- 

We can derive an explicit form for an+i,i in terms of a„.i. 



<ln+l,l 



Un,l 
Un,0 



(2.8) 



where 



oc 
i=0 



(-l)'(2(/ + »))! 
(20! 



p-\-q—l^i 



"n,q 



(2.9) 



B. The high-temperature expansion 

To study the susceptibility not too far from the HT fixed point, we can expand x in terms of We need to expand 
each of the parameters to some power in /? at each level to find %. Since we choose the sc aling factor u so that P is 
eliminated from the recursion, we find that ao,; oc From the form of the recursion, Eq. (|2.9[) , we can see that an,i 
will always have /?' as the leading power in its HT expansion (since p + q > I)- If we want R{k) expanded to order 
/3™""'^ , we will use the truncated recursion formula 

(-1)^(2(^ + 0)! 

Mm^.. = 2^ (A/c)(W)iH2lV ^ l^n,pan.gL_^ , (2.10) 
i=o y I I -y J- p^g^i_^_i 

where the notation [. . .J^^^ means that the expression in brackets should be expanded up to order m„iax in (3- We 
define the coefficients of the expansion of the infinite-volume susceptibility by 
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(2.11) 



m=0 



We define = bm/bm-i, the ratio of two successive coefiicients, and introduce quantities ||T^, called the extrapolated 
ratio (Rm) and the extrapolated slope (S'm), defined by 



and 



where 



Rm = mr,n - (to - l)r„_i , 
Sm = mSm " (to - l)5„i_i , 

^ ^ -to(to - l)(rm - r,n-i) 
mrra - (to - l)r„i_i 

is called the normalized slope. If we use the expansion 

X~(/3c-/3)"^(Ao+Ai(/3c-/3)^ + 
and assume that Aq and Ai are constants, wc find 

5™ = 7 - 1 - isTTO-^ + 0(to-2) 



(2.12) 
(2.13) 

(2.14) 

(2.15) 
(2.16) 



where X is a constant. However, if we calculate this quantity in Z? = 3 with Ising and Landau- Ginsburg measures, 
we find oscillations (Ref. ]l8| , p^ ). For comparison with results described later, Fig. |l| shows these oscillations in the 
extrapolated slope. 
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FIG. 1. The extrapolated slope (S) versus m for the HM with c = 2^/^. 

C. The truncation approximation 

The HT expansion can be calculated to very high order, however, due to a large number of subleading corrections, 
this is a very inefficient way to obtain informat ion a bout the critical behavior. In Ref. ||^, it was found that one 
can obtain much better results by altering Eq. ( 2.101 ) . First, we can retain a much smaller number of terms in the 
sum (originally from 1 to rrimax — than in the expansions [. . .J^^^ . As an example, one can calculate the 1000-th 
HT coefficient of x with 16 digits of accuracy using only 35 terms in the sum. Second, we can simply replace the 
expansions [. . .] by numerical values. This is equivalent to consider the polynomial approximation 



-Rn(fc) 
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(2.17) 
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for some integer Imax at each n-step. The mapping is then confined to an Zmaa;-dimensional space. 

There remains to decide if one should or not truncate to order fc^'^ax after squaring i?„. This makes a difference 
since the exponential of the second derivative has terms with arbitrarily high order derivatives. Numerically, one 
gets better results at intermediate values of Imax by keeping all the terms in R^. In addition, for the calculations 



performed later, the intermediate truncation pa ds th e "structure constants" of the maps (see sec. VII) with about 



fifty percent of zeroes. A closer look at section VII, may convince the reader that not truncating after squaring is 
more natural because we obtain correct (in the sense that they keep their value when Imax is increased) structure 
constants in place of these zeroes. We have thus followed the second possibility where we truncate only once at the 
end of the calculation. With this choice 

(-l)'(2a + z))! ^ 

1=0 ^ ' ' ^ ^ p+q=l+i 

Compared to the HT expansion, the inital truncation to order Imax is accurate up to order /?'""'^ . After one 
iteration, we will miss terms of order but we will also generate some contributions of order /J^'™"^. After n 

iterations we generate some of the terms of order as in superconvergent expansions (such as Newton's method 

to calculate the roots of a polynomial). 



III. A ONE- VARIABLE MODEL 

Before attacking the multivariable expansions of the scaling variables, we would like to illustrate the main ideas 
and study the convergence of series with a simple one variable example which retains the important features: a 
critical temperature, RG flows going from an unstable fixed point to a stable one, and log-periodic oscillations in the 
susceptibility. 

In order to obtain a simple one-variable model, we first consider the Imax 
mapping is then reduced to only one variable. The mapping takes the form: 



<ln+l,l 



(c/2)a„,i - (3cV8)c 



1 truncation using Eq. (2. IS). The 

(3.1) 



1 - (c/2)a„,i + (3c2/16)a2 1 ' 
Expanding the denominator and keeping terms in the mapping only up to order 2 in an+i^i, we obtain 

(c/2)a„,i - (cV8)a2 . 



(3.2) 
(3.3) 



From Eq. (2.6), we can put this recursion directly in terms of the (tnmcation approximated) susceptibility: 

f3 /c\"+i 2 

This approximate equation was successfully used in Ref. Q to model the finite-size effects. If we expand x i'^ 
and define the exptrapolated slope, 5™, as in Eq. ( 2.13D , we see oscillations in Fig. |^ quite similar to those in the HM. 
(Fig. 111). 
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FIG. 2. The extrapolated slope {Sm) versus m for the HT of x calculated from the simplified recursion Eq. (3.J) with 
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Our goal is to obtain the susceptibility (in the one-variable model) as a function of /3 and c/2. Despite the simple 
appearance of the susceptibility recursion, the behavior of the large-n limit is not immediately apparent. Simply 
running the recursion and examining the result yields little, as the expression becomes very complicated in short 
order. Instead, we approximate the susceptibility near to either (3 = or f3 — j3c, then successively build up more 
accurate approximations as we move the temperature away from this point. 

In the following, we use the notation ^ = c/2 and we only consider the case xo — ^ in the Ising model. Our 
analysis will be simplified if we remove the n-dependence of the recursion. For this purpose, we define a new quantity, 
hn, such that 

Xn ^ ^ , (3.4) 
where a is an arbitrary constant. This gives the recursion 



K+i = iK + jS.'^ahl . (3.5) 



The fixed points of this map are 



^*=0'^^^- (3-6) 



We can choose a so that the non-zero fixed point is equal to one: 



(3.7) 



This has the nice effect of making the fixed points independent of ^ and [3. Also, (3 is removed entirely from the map, 
making it only dependent on ^: 

K+i = + (1 - i)hi . (3.8) 

We call this map the "/i-map" . We recover the susceptibility from: 



4(1 - g) K 



Recalling that xo = li we have 



allowing us to write 



Xn = (3.11) 

ho 



for non-zero ho (Hq = means th at /3 — and x„ = 1 for all n) 



If we iterate the recursion Eq. (3.8), we find that for < ft-o < 1, hn as n —^ oo. Correspondingly, for n —^ oo, 
Xn approaches a constant, x, as we saw by using the original recursion. For /iq = 1, /i„ = 1 for all n. For ho > 1, 
h„ oo for n — > oo. This is the same behavior of the susceptibility expected as the temperature (/?) crosses the 
critical value. The adjustable parameter in ho is (3, so we see that ho — I corresponds to f3 — f3c, consequently 

ft ^ ^flp^ . (3.12) 

In Fig . m we show the values of hn for different values of ho- This figure indicates that the shape of the function is 
independent of ho ■ 
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FIG. 3. h-a versus n for C = 0.5. Here /lo = 1 - 10"^ 1 - 10"% and 1 - 10"^% 
To understand the behavior of the ft,-map, we need to investigate the stabihty of the fixed points. We find: 

' dh„A 



dh„ 
and 

dhn+i 



(3.13) 



/i„=0 



dhr, 



= 2-^. (3.14) 



In the foUowing, we only consider < ^ < 1, so the zero fixed point is stable, while the other is unstable. All the 
between zero and the unstable fixed point are attracted toward the stable fixed point. Also note that for all /iq greater 
than the unstable fixed point, ft,„ — > oo. For at the unstable fixed point, hn — 1 for all n. So then Xn = C""- This 
diverges, since ^ < 1. 

Note that the /i-map can be put in the form used in Collet and Eckmann monography pO|] , 

Xn+i = 1 - ^.xl^ , (3.15) 

with 

A* = (e/4)(e-2) , (3.16) 

by using a linear transformation. With their parametrization, the first bifurcation where the fixed point becomes 
unstable and a cycle 2 develops occurs a.i ji — 0.75. This corresponds to ^=3 or -1 which is outside of the region 
where we will study the map in the following (clearly, a cyclic behavior means no thermodynamics limit). 
We can expand the map about the unstable fixed point, hn — 1. 

K+i = 1 + (2 - i){K -!) + (!- OiK - If ■ (3.17) 

The usual notation for the eigenvalue near the critical point is A, so we have A = 2 — ^. If we define 

dn = l-K, (3.18) 

then: 

= Ad„ + (1 - \)dl , (3.19) 

with the starting value do = 1 — (3/ Pc- We call this map the "d-map". 

Note the similarity of the d-map to the original /i-map. We can introduce a duality transformation |lj] between 
the two maps which interchanges /i„ <-> dn and ^ <-> A. If the duality transformation is applied twice, we return to 
the original quantities. For Q < < 1, 'we also have < dp < 1 with small values (approaching from above) in one 
variable corresponding to large values (approaching 1 from below) values in the dual variable. 
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We would like to construct the susceptibility, X: ^ function of do- As a first approximation, we use Eq. (3.19) to 
linearize near the unstable fixed point (do = 0), finding the critical exponent 7 in the process. Beginning with a value 
of do very small (/3 close to /3c), then for a certain number of iterations d„+i ~ Ad„, so that d„ ~ A"do. So as long 
as A"(l — /3//3c) << 1, then /i„ ~ 1 and Xn — C""- If assume that /i„ stays near 1 for some number of iterations, 
and then drops quickly, near some n = n*, to the region where hn oc (thus stabilizing x), then x ~ • 

Defining rt* by A" do = 1, we have 



which gives 



In do 



(3.20) 



(3.21) 



Using the usual notation for the critical exponent, the leading singularity is given by (1 — /3//3c) ^ ■ So we have 



or, equivalently: 



A'^ 



InA ' 



(3.22) 



(3.23) 



We can divide the leading singularity out of x and plot the remainder near to the critical point (Fig. 0). 
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FIG. 4. Oscillations near the critical point. The susceptibility is generated from the recursion Eq. (3.19) with A = 1.9, while 
7 and /3c are calculated from the formulas derived in the text. 

We see periodic oscillations with respect to the variable ln(l — ^) which have period of In A. The linear approximation 
gives no clue as to the origins of the oscillations. To understand this, as well as the higher-order corrections to the 
leading singularity, we need to expand d„ to higher order in do. 



IV. SCALING VARIABLES IN THE ONE- VARIABLE MODEL 

In this section, we show that the idea of scaling variable comes naturally as a way to express the susceptibility of 
the one- variable model in terms of the input parameters. 



A. dn as a function of do 

We can find d„ in terms of do by expanding 
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dn^^qn,idl)- (4.1) 

and plugging into Eq. (2.4). Wc find 



1=1 



i-i 

qn,l = Mn-l,l + (1 - A) ^ qn-l,j qn-l,l-j ■ (4.2) 

Given the initial coefRcients, i.e. qi^i = A, qi,2 = (1 ~ A), and qij = for j > 2, we can find all the coefficients at any 
n. For example, 

,„,, = ^M. ,„) 

Finding forms for the coefficients explicitly in terms of A and n by this direct method quickly becomes difficult. We 
will show that this task can be simplified by introducing the scaling variables. 

As noticed before ^ the transition region from c?„ ~ to dn ~ 1 has a shape which looks independent of do (see 
Fig. H and use Eq. ( 3.18 )). The shape seem to depend only on the value of A. We can prove this in a restricted way 



by imagining a new initial value of "do" which is identical to the current di. The shape of the curve to the right 
would then be identical to the current case. Likewise, we can reverse the map, getting new values for "do" which 
nonetheless generate all the values in the current series do, di, . . .. This inverse map is unique, assuming we confine it 
to only positive values for d„, and has the form 

In this case, the renormalization group is a group in the strict sense, since we have defined a unique inverse. 

What we would like then is a way of parameterizing d„ in terms of a function independent of n. Regardless of 
our actual do, we can extrapolate backwards, as suggested above, to another "do", which is small enough so that the 
linearized method works, where each new "d„" scales approximately as A" times the "do" . This suggests using a new 
parameter, which we call j/„, that scales exactly like this, so that if d„ corresponds to then d„+i corresponds to 
Un+i = Aj/„. This is our first covariant quantity, so-called since it scales exactly the same, regardless of the value of n. 

The idea of using variables with simple transformation properties has a long history, for instance the angle-action 
variables in classical mechanics and the normal form of differential equations appearing in Poincare's dissertation pl| . 
For continuous RG transformations, Wegner introduced the notion of "scaling field" . We discuss here its analog 
for discrete RG transformations. 



Let us define a function d such that d(y„) = d„. From the explicit form of the recursion formula Eq. {2A) this 
requirement implies 

d{Xyn) = Ad(y„) + (1 - X)d^{yn) ■ (4.5) 

If we let 



d(2/„) = ^s,yi , (4.6) 

we find that 



1=1 



1-x 



si 



A' - A . ^ 



J2 si-j . (4.7) 



The first coefficient is undetermined. We let si = 1, so that for small j/„, we have d„ — Un- The first few coefficients 
give 

d{yn) = yn- jvl + x^^x + lf "" ^ ■ ■ ■ ^^'^^ 



In Fig. H, we plot d(y„) for several values of A. 
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FIG. 5. d{y„) for various values of A. The top and bottom curves represent the A — > 2 and A — > 1 hmits, respectively. The 
dash-dot-dot line is A = 1.9, the dotted line is A = 1.6, and the dash-dot line is A = 1.3. 

In the limit where A ^ 2, we can calculate explicitly d(y„). This corresponds to ^ = 0, so that 

hn+i = hi . (4.9) 

This is satisfied if hn — e°^", where a is arbitrary. So we have c?„ = 1 — e°^". To satisfy yn — dn for small d„, we 
need a = — 1. Finally, we have 

d(y„) = 1 - e-y- . (4.10) 

The expression of c?„ as d{\"'yQ) allows us to interpolate between integral values of n. In Fig. ^ we show the curve 
for d{yn) superimposed on points generated from the c?-map for A = 1.5. 
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n 

FIG. 6. dijjn) plotted against n. The data points (circles) are generated directly from the d-map recursion with A = 1.5 and 
do = 10"*. 

We can invert the series for d{yn) in terms of y„, defining a new function y, so that y{dn) = yn- Since we know the 
coefficients of the d{yn) series, we can invert to get y in terms of d„. More directly, we can use yn+i = Aj/„, or 

y{\dn + (1 - \)dl) = \y{dn) , (4.11) 

along with the constraint y{dn) — dn for small values of to construct a series solution. If we let 

OQ 

y{dn)^J2^idl , (4.12) 

1=1 
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we find that 

int{l/2) _ 

with ti = 1. The notation int{l/2) means to use the integer part of 1/2. The first few terms give: 

y{dn) = dn + jdl + j^^—^di + ... (4.14) 

Taking our A — > 2 example, we can immediately invert our result for d{yn), to find 

y{dn) = -H'^-dn) . (4.15) 
With the t/-function, we can now reconstruct d„ in terms of do: 

dn = y-\y'y{do)) . (4.16) 



B. /i„ as a function of ho 

Because of the duality between our two maps, wc can easily reproduce all of the above results of the d-map for the 
/i-map. Just as with dn, we can parameterize hn in terms of a covariant variable, y„, dual to y„. Similarly to before, 
we define h{yn) = hn, with yn = ^""yo- We can also define the inverse function for h, y, where 

mhn + il-Ohl)=mhn) ■ (4.17) 

We can immediately carry over all of the results we found for d{yn) and y{dn) to /i(y„) and y{hn), replacing A with 
^. For example, 

1 2 

KVu) =yn- jfn + ^2(^ + l) ^n + " " " ' (4-18) 



while 



y{K) = K + jhl + i^^^^n + • • • (4-19) 

As with dn, we can find /i„ to any order in Hq. The first few terms in the expansion give: 

hn = Cho+ ^ hl + ... (4.20) 

The new covariant variable behaves similarly to y„. However, since ^ < 1, we see that as n ^ oo, y„ 0, while 
as n — * —00, yn oo. As with the d-map, y„ goes to zero with while blowing up as hn approaches 1. Thus, for 
various ^, h versus y„ and y versus hn will have at least qualitatively similar shapes to the corresponding d-map plots. 

V. LARGE ORDER BEHAVIOR FOR THE ONE- VARIABLE MODEL 

In this section we give empirical results concerning the large order behavior of y{d), y{h) and their inverses. From 
this, we infer the domain of convergence of the expansions. We would like to know if the series converge for the ranges 
of interest. The allowed ranges are to infinity for and y„, and to 1 for d„ and /i„. 
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A. d{y) 



For the d-map, we first e xam ine tlie extreme allowed values for A, 1 and 2. For A = 2 = 0), we already found 
that d{y) = 1 — e ~^. Eq. (4/7) converges in the whole complex plane. At the other end of our range, for A ^ 1 
^ 1), Eq. dO] ) becomes 



-1 



si 



(5.1) 



so that 



d{y) =y-y^ + y''-y^ + ... (5.2) 

This converges only for < y < 1, to d{y) = y/{l + y). 

We use the ratio test to determine convergence for our series. The series will converge in the entire complex plane if 
|si_i/s/| increases without bound as I oo. For all tested values of 1 < A < 2, we found that we could fit a straight 
line to plots of In |s/_i/s/| versus \n{l) for large enough I (Fig. R). 




t>.5t) 5.6 b.6b 5.7 



ln(l) 

FIG. 7. Linear fit to a logarithmic plot of the ratio of coefficients against the order of the bottom coefficient. Here, A = 1.5. 
The slope of the line is 0.584. 



Thus for large Z, the coefficients follow the rule 

\^\^Cl\ (5.3) 

Sl+l 

where we find that C is always of order 1 (in fact, 0.9 < C < 1) and < fc < 1. The value of k can be determined by 
a linear fit. It seems clear that with 1 < A < 2, the series for d{y) converges for all < y„ < oo. In agreement with 
the above exact series, for A 1, C — > 1 and A: 0, while for A ^ 2, C ^ 1 and k ^ 1. The intermediate values are 
given in Fig. ^. 
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FIG. 8. The exponent k defined in Eq. (5^) as a function of A 



The form that the ratio of successive coefficients takes also gives us an idea of the number of coefficients we will 
need to approximate d{y) for a given value of y. Assuming that a particular term in the series becomes significant 
when it is the same size as the previous term, we find that the maximum order that we need in our approximation 
goes like 



l^y 



l/k 



(5.4) 



Thus for A near 2, the number of terms needed will grow approximately linearly with y. On the other hand, for A 
near 1, where fc — *■ 0, the number of terms will increase much more quickly with y. 



B. h(y) 

We now turn to the dual function, h. Above, we found the series in d(y) for A ^ 1. From this, we can immediately 
write down h{y) when ^ ^ 1. The series is identical in form, with A replaced by ^, so that 

Ky)^y-f + f + ... , (5.5) 

which converges for < y < 1 to h(jj) — y/{l + y). For the general case, we again examine the ratios of successive 
coefficients, |s//si+i|. We find the ratios flatten to constant values, for large enough I , indicating finite radii of 
convergence. The radii get smaller and vanish as ^ approaches zero as shown in Fig. 0. 
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FIG. 9. Radius of convergence of h{y) as a function of A. 

Note that for any value of ^ tried, we found good evidence that limn—,oaSn = —1.5, indicating a {y — ycY^^ behavior. 
This is consistent with the existence of a quadratic minimum for the inverse function discussed below. 



C. y{d)) 

We already found above that for A = 2, y(d) = — ln(l — d). This series converges for < d < 1, which is the 
required range. (Since y must diverge as d ^ 1, it is not surprising that our series breaks down at this point.) For 
A ^ 1, we have, inverting the c?-function found above, that y{d) = d/{l — d) which also converges for < d < 1. We 
find empirically from the analysis of ratios that for all 1 < A < 2, converges in the region < c? < 1. 

The analysis of the extrapolated slope for various A gives convincing evidence that the main singularity has the 
form 

yid) ^ d-^h . (5.6) 
This can be seen with short series when A is close to one and requires larger and larger series as A gets close to 2. 



D. m) 

The y{h) series behaves similarly to the y{d), having a radius of convergence of 1 for all ^. In this case, the ratio 
of coefficients ti/ti+i approaches 1 smoothly from below as / increases, for larger ^ Around ^ = 0.4, small oscillations 
begin to appear in the curve. As ^ becomes smaller, these oscillations grow and become quite large for ^ smaller than 
0.1 (Fig. |l^). However, in all cases studied, the ratio eventually approaches 1, for large enough 
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FIG. 10. Ratios of coefficients for y{hn) with ^ = 0.1. 

A detailed study shows that if we continuate y{h) for negative values of h using the series expansion, the function 
develops a quadratic minimum when h —1. The absolute value of this minimum found in all examples studied is 
compatible with the radius of convergence of the inverse function and the square root singularity of the h{y) discussed 
above. 

As one can guess by looking at Fig. |2| the analysis of the extrapolated slope is intricated. However, if we calculate 
enough terms and if A is not too close to 2, one can get approximate results which are consistent with a main singularity 
of the form 



y{h) ^ h- 



(5.7) 



For instance, just by looking at the asymptotic behavior of Fig. one can see that 7 ~ 1.4677, as expected, with 
errors of the order 10~^. 



VI. THE SUSCEPTIBILITY OF THE ONE- VARIABLE MODEL 



Now that we can calculate /i„ and (i„, we can construct the (infinite- volume limit) susceptibility. If we focus on the 
HT fixed point and the /i-map, we get the HT expansion in terms of ho = (3 / 13c- On the other hand, to understand 
the susceptibility's behavior near to the fixed point, we need to expand in terms of do = 1 — (3/ 13c- Finally, we consider 
the possibility of combining the two expansions in the crossover region. 



A. The HT expansion of the susceptibiHty 



We could find the HT expansion from the finite-volume susceptibility found in the previous section, by taking the 
large-n limit. Since < ^ < 1, all of the ^" terms will drop out. However, we can also construct the large-n limit 
directly. Recalling that Xn — (^n//io)C~", and using y„ — ^"yo, we find that 

_ Kyg _ h{yn)yo . , 

hoyn hoy-a 



The infinite- volume limit becomes: 



Using the y-expansion for h, we find 



X = lim Xn = 7- 1™ — : — . (6.2) 

- ho n^oo y„ 



n — ^00 



Vn C 

As 71 ^ cxD, jjn 0, SO the infinite volume limit of the susceptibility is: 



^^(^"^-1-% + ... (6.3) 
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x=yj^ = yM, (6.4) 

ho ho 

Using the /i-expansion for y, we get the HT expansion for the susceptibility: 

X=l + ^ho + ^^^—^hl + ... (6.5) 

Using Xn — (^n/^o)f~", we find the HT expansion for the finite-volume susceptibihty (recall that ho = (3/ Pc)- 

1 - f " 

X„ = 1 + -^ho + . . . (6.6) 

We can find the error in the finite- volume susceptibility as an approximation to For large n (small y„), we get 
Xn — x(l - Vn/O = x(l ^ r'""^yo)- The relative error is then: 

f-iyo = r^'/^oX • (6.7) 



X 

This behavior is observed H| in the approach to infinite volume for the HM. 



B. Expansion near the critical point 

In addition of using the HT expansion (expressing x i^i terms of ho), we would also like to have an expansion in 
terms of the dual variable do. Multiplying the numerator and the denominator of Eq. (^^) by ?/q , we obtain 



" (1 - dMd„)i 

with 

e = yovl ■ (6.9) 

But since — , we have also 

e = vuvl , (6.10) 

for any n. In other words, is RG invariant. Due to the discrete nature of our RG transformation, Q is not exactly 
a constant. If expressed in terms of In(yo), 6 is a periodic function of period InA. However, for A not too close to 2, 
the non-zero Fourier modes are very small um . 



C. Numerical evidence for overlapping domains of convergence 

As we have seen above, ynUZ is the same for every n. We can thus pick n such that we are just in the crossover region 
and both expansions have a reasonable chance to be accurate. In order to test the accuracy of the approximations 
Vappid) (series expansion up to a certain order) for various n, we have prepared an empirical sequences of dn starting 
with do = 10~^. We have then tested the scaling properties by calculating 

Dn = I [yappidn)/iyo\n] - M , (6.11) 

where yo was calculated with 16 digits of accuracy. Optimal approximations are those for which £)„ ~ 10^^^. For such 
approximation, scaling is as good as it can possibly be. Indeed, due to the peculiar way numerical errors propagate 



^2| , one does not reach exactly the expected level 10"^^ (more about this question in section IX). We can define a 
similar dual quantity by replacing d by ft. and y by y. In this case, yo is estimated with the same accuracy as yo by 
stabilizing yihn)/^^, for large enough n. 

We have performed this calculation for A = 1.1, 1.5 and 1.9. The conclusions in the three cases are identical. For n 
large enough, the D„ of y starts increasing from 10~^^ until it saturates around 1. By increasing the number of terms 
in the expansion, we can increase the value of n for which we start losing accuracy. Similarly, for n low enough, the 
Dn of y starts increasing etc... We want to know if it is possible to calculate enough terms in each expansion to get 
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scaling with some desired accuracy for both functions. The answer to this question is affirmative according to Fig. 
|ll|for A = 1.5. 



1 

0.01 

1E-4 

1E-6 

1E-8 
Q 1E-10 
1E-12 
1E-14 
1E-16 
1E-18 

30 35 40 45 50 55 60 
Iteration Number n 

FIG. 11. Departure from scaling D„ defined in the text, for y (curves reaching 1 to the right) and y (curves reaching 1 to 
the left). In each cases, we have considered approximations of order 10 (dot-dashes), 30 (dots), 50 (dashes) and 70 (solid line). 
The value of A is 1.5. 

One sees, for instance, that with 10 terms in each series, we have scaling with about 1 part in 1000 near n = 45 for 
both expansions. The situation can be improved. For 70-70 expansions, an optimal accuracy is reached from n — 44 
to 46. For the other values of A quoted above, similar conclusions are reached, the only difference being the optimal 
values of n. 



D. Analysis of the RG invariant 

Another evidence for overlapping convergence is that we can calculate the RG invariant O for a certain range of 
yn- To evaluate O, we use the series expansions for y and d, cutting each off at some order: 

ni m 

y{l-d{y))^Y.*^{l-Y,s,y^y , (6.12) 

i=i j=i 

where si and ti are the Ith coefhcients in the d and y series, respectively. The d series is more accurate the smaller is 
the value of y, while the y series is more accurate for d{y) close to 1, which means large values of y. Thus for y very 
large, we need many terms in the d series, while we need less terms in the y series. On the other hand, when y is very 
small, we need a large value for rh and a relatively small m. We have found that, given a fixed value of m -\- m, the 
most accurate values for 8 are obtained when m ~ m. In Fig. |l^, we show 8 calculated by keeping 50 terms each in 
the expansions for y and y. The result is plotted against ln(?/). We used ^ = O.I, which makes the oscillations much 
larger than, for example, near to ^ = 2~^/^. Near the fixed points, we need more terms in the appropriate series to 
get accurate results. 
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In(yo) 

FIG. 12. The invariant function O, calculated at ^ = 0.1, and plotted against ln(i/o) 



We can study the oscillation we see in Q by looking at its Fourier expansion. Since O is periodic in Inyoi we can 
express 



0(2/0) = VoVil ~ d{yo)) = ^ap 



^ipu In yo 



(6.13) 



where u; = 2tt/ In A. The coefficients are given as 



1 

i^Ta 



(6.14) 



for some ya of our choice. 

As an example, we calculate aq for ^ — 0.1 where the oscillations are substantial. The choice of the interval can 
be infered from Fig. If we had infinite series, the function would be exactly periodic. For finite series, we see 
that ya cannot be too large or too small. For intermediate values, we obtain oq ~ 6.06676. In Fig. |l^, we show this 
constant term subtracted from the Q we evaluated above. What remains is oscillations about 0. We also calculated 
fli (using ya — X so that we keep information about the phase constant as well as the amplitude). We construct the 
first oscillating mode from this, and in Fig. ( p^ we show the result of subtracting out this in addition to the constant 
term. We are left with even smaller oscillations, from the 02 term. 




in(yo) 

FIG. 13. Same as previous plot, with the constant term from the Fourier expansion subtracted out. 
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FIG. 14. Same as previous plot, with first oscillating term from the Fourier expansion subtracted out. 

In Fig. |l^. one notices that when In(yo) is of the order of 1.7, the breakdown of the approximation starts with small 
numerical noise before the complete breakdown near ln(?/o) ~ 2.2 is observed (see Fig. n2h. 



VII. SCALING VARIABLES IN THE HIERARCHICAL MODEL 



In the one- variable model, we learned that we could express the basic mapping parameter, /i„, the distance from the 
HT fixed point, in terms of either of two new variables, ?/„ or y„. The first of these variables transforms multiplicatively 
with A, the eigenvalue near the unstable fixed point, while the second transforms multiplicatively with ^, the eigenvalue 
near the HT fixed point. 

We would like to extend these methods to Dyson's HM. We will express the parameters of the model in terms of 
two new sets of variables yn,i, ■ ■ ■ :yn,imax ^-'^d yn,i, ■ ■ ■ ^ynj^a^j where Imax is equal to the number of parameters we 
keep in the mapping. These are the scaling variables. 

Near to the unstable fixed point we find a set of eigenvectors and eigenvalues which control the RG flows in the 
linearized region. We use the notation for the ith eigenvalue, ordered in size from greatest to least. For each of these 
eigenvalues, we introduce a variable yn^i which transforms under the RG transformation as yn+i,i = Kyn.i- We denote 
Xi for the eigenvalues in the linear HT region, similarly defining variables which transform as yn+i,i = \yn,i- Then 
we can expand each of the parameters in terms of each set of scaling variables. Our first task will be to use special 
coordinates where the problem will be solved at the linear level. To use a relativistic analogy, the transformations 
leading to such a system of coordinates (the "linear scaling variables"), can be interpreted as a "translation" or shift 
and a "rotation" or more generally a linear transformation. 



A. Transformation properties under translation 



We previously considered truncation in the an,i parameters. Unfortunately, this approximation becomes less ac- 
curate if we are interested in the critical region. We would like to investigate a truncation that works near to the 
nontrivial fixed point. We re-express the basic recursion in terms of a new set of parameters, which give the distance 
from the nontrivial fixed point. For notational convenience, we first rewrite the unnormalized recursion given in Eq. 



( 2.18 ), using the "structure constants": 
with 

f (r/A\t^+'^ {-l/2r+ — {2{^c+^)y. for/, + ,y>n- 

= < ^ ' > {fJ^+<y-a)\(2a)\ ' tor + ^ a 2) 

'^1 , otherwise . 

These zeroes can be understood as a "selection rule" associated with the fact that an,i is of order as explained in 
section ||. If we follow the truncation procedure explained in section ||, the indices simply run over a finite number 
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of values. We use "relativistic" notations. Repeated indices mean summation. The greek indices indices jj, and v go 
from to Imax, while latin indices i, j go from 1 to Imax- Obviously, V^"^ is symmetric in fi and v. By construction, 
o-n,o — 1 and we can write the normalized recursion in the form: 

with 

7w; = 2r°' . (7.4) 

Let assume that we know a fixed point of the recursion formula an,i = a* . We then introduce intermediate variables 
Qnj which are zero at the fixed point: 

an,i = a'l + gn,i ■ (7.5) 



Plugging this relation into Eq. (7.3), using the fixed point properties, subtracting the fixed point and reducing to the 



same denominator, we obtain a recursion formula for the gn i having the same form as Eq. (7.3) with the substitutions 



r;'-^;?7(r;-''-(i-M«;n^) , (7.6) 





and 

- AiM], + 2r^^a* + (1 - 6,,){Ml + 2Tl^a;)a; , (7.7) 

where Uq = TQ'^a*^a*. 

These equations express the transformation properties of the structure constants under a translation (shift) of the 
coordinates. In particular if a* = (as for the HT fixed point), one can see that since Fqq = 1, the transformation 
reduces to the identity. We will now study the transformation properties under changes of coordinates which will 
diagonalize the linear RG transformation matrix . 



B. Diagonalization of the linear RG 

The diagonalization of the linear Al^ near the HT fixed point is quite simple because it is of the upper triangular 



form. The eigenvalues are just the diagonal terms. From Eq. (7.2), one sees that Ith diagonal term is given as 



A(,) = 2(c/4)'- . (7.8) 

This spectrum was obtained in Ref. |p3| with a different method. This means that an ^ x ^ truncated version of 
the matrix will have eigenvalues which are identical to the first I eigenvalues of the full HM map. Furthermore, the 
Zth eigenvector will contain only I non-zero entries. This means that if we truncate to a;^^^ , we are simultaneously 
truncating to the parameter space to a subspace spanned by the first Imax eigenvectors. This "stability" is due to the 
special relationship existing between the HT expansion and the polynomial truncation explained in section ||. 

On the other hand, the linearized map near the nontrivial fixed point does not have these simple properties. The 
eigenvalues may be determined numerically from a truncated version of the linearized matrix. With a large enough 
truncation, we can find a certain number of the eigenvalues to any desired precision. We find that there is one 
eigenvalue larger than 1, with all the rest less than 1. Though we do not have a simple formula for these eigenvalues, 
we know that they do shrink in size quickly, similarly to the eigenvalues of the HT fixed point. For instance the 
numerical values for c = 2^/-^ are Ai = 1.42717 . . ., A2 = 0.85941 .... A more complete list is given in Ref. | pl| . 

Since M is not a symmetric matrix the left and right eigenvector are distinct. 

Mir, = A(,)Vr (7.9) 

and 

4MI = A(,)0; . (7.10) 
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The notation (r) means that there is no sum on r. Since all the eigenvalues are distinct, one can normalize the 
eigenvectors in such a way that 



4>i^: =51 . (7.11) 

Similarly a completeness relation (decomposition of the identity) can be obtained by summing over r: 

= 5] . (7.12) 

In numerical calculations, the relations of orthogonality and completeness provide a reliable way to check the accuracy 
of the calculations. 

We can thus diagonalize the linearized RG: 

4,[M\i,i = X^r)S:' . (7.13) 



If we start with the transformation expressed in the general form Eq. (7.3) but with a replaced by g to indicate that 



the fixed point is at the origin, as in the previous subsection, we can define a new set of coordinates dn,i, such that 

gn,l = '>Pldn,r , (7-14) 

or, equivalently, 

dn,r = (t>ign,l ■ (7-15) 

With these notations, all the structure constants transform according to the familiar rules of tensor analysis. Since (p 
and V' are inverse of each other, we ca n thin k of lower indices as covariant and upper indices as contravariant. As an 



example of transformation under Eq. ( 7.14 ) , we have 

n^^^vyr^.r,, . (7.16) 

The other structure constants transform according to the same rules. 

The eigenvectors are not determined uniquely even after normalizing. If ipl and are components of normalized 
right and left eigenvectors, then ipl = a(^r)^i and = 4)l./a(^r)i where a(r) is a constant, will work equally well. In the 
following, we will fix this ambiguity by requiring that the "other" fixed point be at (1, 1, . . .) in the new coordinates. 

C. The Canonical Form 

In summary, we can choose a system of coordinate di where the unstable fixed point will be at the origin of the 
coordinate and the HT fixed point at (1, 1, . . .), and such that the linear RG transformation is diagonal. In this system 
of coordinates, the RG transformation reads 



X(r)dn.7^ ^" dji^pdji q 

1 + APdn,p + Ag-^ 

dn.pdn,q 



a„+l,r — , ^„ J , Apq J J J yi-^l) 



where the new structure constant are calculated from the original ones according to the transformation laws discussed 
in the previous two subsections. 

Similarly, we can introduce new coordinates hi, so that 

and = i'\hn,r , (7-18) 

and take the eigenvectors so that the unstable fixed point is at (1, 1, . . .) in the new system of coordinates. Note that 



the form of the eigenvectors guarantees that /i„j is of order /3'. This can be seen by inverting Eq. ( 7.1§| ) using the 
matrix of left eigenvalues. Due to the upper-diagonal form, the second left eigenvector has its first entry equal to 
zero, the third its first two entries etc... . In this new system, the RG transformation reads 

hn+i.r = r— Apqu T — ■ -ly) 
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D. Expansions in terms of the scaling variables 



We would like to express the evolution of the canonical coordinates in terms of functions of the scaling variables: 

dn,r= ^r,i,i^...yn\iyn,2- ■ ■ ^ C^-^O) 



where the sums over the z's run from to infinity in each variable and yn^i — A";-)yo,;- In practice, we expand in terms 
of some finite number of the scaling variables, less than or equal to the number of d parameters we have truncated 
to. Using the notation i = (ii, Z2 • • •) and the product symbol, we may rewrite the expansion as 

du,r=J2'r,iY[yl:^ (7.21) 



Using the transformation law for the scaling variables, we have 

dn+l,r = gr,i ]~[(^(m)ym)''" ■ (7.22) 



Each constant term, Sr,o,o,...-: is zero, as the scaling variables vanish at the fixed point. From Eq. (7.17), we see that 
all but one of the linear terms are zero for each value of r. The remaining term is the one proportional to the rth 
scaling variable. We take these coefficients to be 1, so that the dn^r — yn,r for small yn.r- For the higher-order terms, 
we obtain the recursion 

_ J2j+ki=ii'^r'^^P,3^qM ~ -^rj rim '^(m)^^'^P,k) — Sj+k+l=i '^''J Urn ^^{m)^o'^ ^PM^<1^ , ^ 

(rim -^(m) ~ j 

The calculation can be organized in such way that the r. h. s. of the equations are already known. This will be the 
case for instance if we proceed order by order in X^g*?' the degree of non-linearity. As one can see, this expansion 
may suffer from "small denominator problems" . This issue will be discussed elsewhere |2j] . 

We can likewise expand each h^^r in term s of scaling variables yn,\-,yn,2-, ■ ■ ■■ The derived recursions are identical 



in form to those derived above. From Eq. (7.8), one sees that the denominator will vanish for some equations and 
unless the numerator is also zero, the expansion is ill-defined. Again this question will be discussed elsewhere | |2^ . In 
the following, we will only use the expansion for the d-functions. 



E. Expansion of the scaling variables 

One can likewise find expansions of the scaling variables in terms of the canonical coordinates, by setting 

yn,r=J2^r,l]l'^n:m. (7-24) 
i m 

and requiring that when d„ is replaced by d„+i, the function is multiplied by A(r). Since d„+i has a denominator, it 
needs to be expanded for instance in increasing order of non-linearity. The same considerations applies for y. Note 
that small denominators may also be present in these calculations. However, the calculation of yi and yi is free of 
such a problem since the largest eigenvalue cannot be written as product of lower eigenvalues smaller than 1. 

It is in principle simple to obtain these expansions order by order in the degree of non-linearity, however there 
exists some practical limitations. For instance, if we want to calculate all the non-linear terms of order 10 in any of 
the variables, with lmax—3Q, we need to calculate and store 30^*^ ^ 10^^ numbers. In the next section, we show that 
one can organize these calculations in a way which allows accurate answers for the susceptibility. 



F. The susceptibility 



From Eqs. (2.6) and (7.18|), we obtain 



-(2//3)V[/i„,.(2/c)^ 



(7.25) 



22 



For n large enough, the hnear behavior applies and the the hn^r get multiplied by 2(c/4)'" at each iteration. In the 
large n limit, only the r = 1 term survives and consequently, 

X = -(2//3)i^ilim„_o,/i„,i(2/c)" . (7.26) 

Using the same method as in the one-variable model, we can in the limit replace hn^i by and obtain 

X - -(2//3)^iyo,i . (7.27) 

This allows us to write 

X = -{2/(3)^ley^ , (7.28) 
where 7 = — In(Ai)/ In(Ai) and O the RG invariant 

© = yoAVo,! ^ yn,iy,1,i • (7-29) 

For reference, in the case c = 2^^'^, ipl ~ —0.564. 

One can calculate the subtracted four-point function following the same procedure, namely expressing the an,i in 
terms of the hn,i- However, > A2 and one needs to go beyond the linear expansion to calculate these quantities. 



VIII. APPROXIMATIONS AND NUMERICAL IMPLEMENTATION 



In this section, we show that it is possible to design approximations such that one can calculate the susceptibility 
using Eq. (7.29). For this purpose we have calculated an empirical series of with c = 2^/'^, (3 = [3c — 10~* and an 
initial Ising measure. Detail relevant for this calculations can be found in Refs. |^,|ll|,^. The calculations have been 
performed with Imax ~ 30, a value for which at the (3 considered, the errors due to the truncation are of the same 
order as those due to the numerical errors. These errors are small enough to allow a determination of the susceptibility 
with seven significant digits if we use double precision. 

The empirical flow proceeds in two steps. First, the flow goes from the initial mesure to close to the unstable fixed 
point. Second, the flow goes from close to the unstable fixed point to the HT fixed point. The first step depends 
on the choice of the initial measure and will not be discussed in full detail. Our main goal will be to show that it is 
possible to construct nonlinear expansions which allow to describe accurately the second step. We now discuss the 
flow chronologically. In order to keep track of the chronological sequence, we postponed the technical discussion of 
the convergence of the series to section [X. In this section, the results are simply quoted. 

Our choice of (3 means that we start near the stable manifold. After about 25 iterations, we start approaching 
the unstable fixed point and the linear behavior dn+i,i — A;(i„ becomes a good approximation. During the next 20 
iterations, the irrelevant variables die off at the linear rate and at the same time we move away from the fixed point 
along the unstable direction, also at the linear rate. At n = 47 we are in good approximation on the unstable manifold 
and dn^2 becomes proportional to d\^. In other words, the non-linear terms are taking over. At this point, we can 
approximate the dn,i as functions of yi only: dn,i — (i;(A"2/i, 0, 0, . . .). This approximation is consistent in the sense 
that if ?/2 = at n = then it is also the case for positive n. 



We have calculated (ii(?/i,0,0, . . .) up to order 80 in j/i using Eq. (7.23). We have then inverted di to obtain j/i. 
Given the empirical dn^i we then calculated the approximate yn,i and then used the other functions di{yi) (with / > 2) 
to "predict" Gomparison with the actual numbers were good in a restricted range. For n — 49, the relative errors 
were less than a percent. They kept decreasing to less than one part in 10,000 for n = 54 and then increased again. 
It will be shown in section IX that this corresponds to the fact that when yi becomes too large (a value of 3.7 first 



exceeded at n = 57), the series expansion of di diverges, unlike the one- variable model for which d{y) is analytical. 
It will also be shown later that the quality of the approximation between n = 45 and n = 55 can be improved by by 
treating 7/2, 2/3, •■ • perturbatively. 

Near n = 55, the presence of the HT fixed point starts dominating the flow but we are still far away from the linear 
regime. To take into account the non-linear effects, we have calculated yi(h) up to order 11 in /?. The reason for this 
modest order is that this is a multi-variable expansion. Recalling the discussion about the HT expansion in section 
we can count the number of terms at each order in /?. At order two, we have /if and /12, but since the linear 
transformation is diagonalized, /12 will only appear in 2/2 (h). At order three, we have hf, /12/11 and not /13. It is easy 
to see that at order Z, one has p{l) — 1 terms, where p(l) is the number of partitions of I. It has been known from the 
work of Hardy and Ramanujan that 
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(8.1) 



It seems thus difficult to get very high order in this expansion. In order to fix the ideas, there are 41 terms at order 
11, 489 at order 20 and 133|8,649 at order 80. 

As explained in section IX, the expansion up to order 11 has a sufficient accuracy to take over at n = 55. It also 
provides optimal (given our use of double precision) results for n > 60. As n increases beyond 60, one can see the 
effect of each order disappear one after the other. Finally, the linear behavior becomes optimal near n = 130. This 
concludes our chronological discussion. The main point made was that there exists a region in the crossover where 
both expansions are valid. We now proceed to justify empirically the claims made above. 



IX. LARGE-ORDER BEHAVIOR FOR THE HM 



A. yi{di) 

The behavior of j/i(rfi) obtained by the procedure described above, has been studied following the same methods 
as for the one-variable model. In order to provide a comparison, we have calculated the same quantities for the 
one- variable model with A = 1.2573. In the following, we call this model the "simplified model" (SM). For this special 
value, the critical exponents 7 of the two models coincide with five significant digits. We have good evidence that 



both series have a radius of convergence 1 as indicated by the extrapolated ratio defined in Eq. ( 2.12 ) reaching 1 at 
an expected rate (Fig. nF 




Order m 



FIG. 15. The absolute value of the diflFerence between the extrapolated ratio and 1 for the HM (empty boxes) and the SM 
(full circles), as a function of the order. 

Similarly, their extrapolated slopes seem to converge to the same value I/7 — 1 ~ —0.23026 ... as shown in Fig. 111. 
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FIG. 16. The extrapolated slope Sm for the HM (empty boxes) and the SM (full circles) as a function of the order . 
In conclusion, the function y{d) for the SM is a reasonably good model to guess the asymptotic behavior of yi{di] 

B. di(2/i) 

The situation is different for the inverse function di (yi ) . A first look at the differences is given in Fig. 
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FIG. 17. Logarithm of the absolute value of the coefficients of the expansion of di(yi) divided by the order, for the HM 
(empty boxes) and the SM (circles). 

The quantity plotted in this figure will be used to discriminate between a finite and an infinite radius of convergence. 
If \bm\ ~ as for a radius of convergence R, then we have ln(|6m|)/TO ^ — In(i?) + A/m for some constant A. On 
the other hand, if |5„j| ^ (m!)~" as for an infinite radius of convergence, then we have ln(|6m|)/"^ ^ — Q;(ln(m) — m). 
In the following, we will compare fits of the form Ai + A2/m and Bi\n{'m) + B-z- 

We first consider the case of SM, where according to the our study in section we should have an infinite radius 
of convergence. This possibility is highly favored as shown in Fig. 
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One sees clearly that the solid line is a much better fit. The chi-square for the solid line fit is 200 times smaller. In 
addition Bi ~ —B2 as expected. In conclusion, this analysis comfirms the ratio analysis done previously and favors 
strongly the inifite radius of convergence possibility. 

The analysis for the HM is more delicate. One observes periodic "dips" in Fig. ^ which make the ratio analysis 
almost impossible. We have thus only considered, the "upper envelope" by removing the dips from Fig. The fit 
represents an upper bound rather than the actual coefficients. The fits of the upper envelope are shown in Fig. |l^ 
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FIG. 19. Comparison of fits of the form Ai + A^jm (solid line)and B\ti\(rn) + B2 (dots) with selected points of the data in 
Fig. |l| for the HM (boxes). 

The possibility of a finite radius of convergence is slightly favored, the chi-square being 0.4 of the one for the other 
possibility. Also, the second fit does not have the B\ ~ —B2 property. From A\ ~ —1.32, we estimate that the radius 
of convergence is about 3.7. 
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C. y^{h) 



As explained in section VIII, one can calculate yi(h) using an expansion in As in section we will use an 
empirical series a„_;, calculate the corresponding hn,i an plug them in the scaling variables. This empirical series was 
calcu lated with an initial Ising measure and f3 = Pc — 10~* (see Ref. Again we define a quantity Dn as in Eq. 

( |6.1lD which is verysmall when we have good scaling and increases when the approximation breaks down. The results 
are shown in Fig. for successive orders in (3. 




FIG. 20. The quantity D„ defined in Eq. (3.11) for expansions of yi in (3, at order 1 (solid line), 2 (dashes), 3 (dots), etc... 
for each iteration n. 

The solid line on the right is the linear approximation. It becomes optimal near n = 130. The next line (dashes) 
is the second order in /3 expansion. It becomes optimal near n — 90. Each next order gets closer and closer to be 
optimal near n = 60. The last curve on the left is the order 9 approximation. It is hard to resolve the next two 
approximations on this graph. 

The asymptotic value is stabilized with 16 digits and one may wonder why we get only scaling with 14 or 15 digits 
in Fi g. |20| . The reason is that we use empirical data and that numerical errors can add coherently as explained in 
Ref. |22[1 . This can be seen directly by considering the difference between two successive values of Dn as shown in 
Fig. |2lf 
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FIG. 21. Difference between two successive D„ in 10"^^ units. 

One sees that the numerical errors at each step tend to be negative more often than positive, and consequently there 
is a small "drift" which affects the last digits. 



D. overlapping domains of convergence 



We can now look at the D„ defined as in Eq. (6.11) for yi and yi together in Fig. Ea 




FIG. 22. Values of D„ for yi up to order 5 (filled circle) and 11 (filled boxes), and for yi up to order 10 in di (empty circles) 
and with first order corrections in y2 (dots), and up to order 80 in di (empty stars) and with first order corrections in 1/2 (solid 
line). 



If we use an expansion of order 5 in (3 for y and of order 10 in di for j/i, we can get scaling within a few percent for 
both variables at n — 54. We can go below 1 part in 1000, with an expansion of order 11 in (3 and order 80 in di. At 
this point, the main problem is that the effects of the subleading correction makes the scaling properties worse when 
n < 57 and n decreases. One can improve the scaling properties by taking the effects of y2 into account. A detailed 
study shows that one can estimate the subleading effects between n = 40 and n = 45. One finds that 
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^i^^ ~ 7.2778 X 10-9 - 4.5 X 10-1° X (9.1) 
A" 

It is thus possible to get a function scaling better by subtracting these correction. This improve the scaling properties 
by almost one order of magnitude near n = 54 and by almost two order of magnitude near n = 45. It is likely that 
our approximation of having a single scaling variable can be corrected order by order in ?/2, 1/3 etc... . 



X. CONCLUSIONS 



We have shown in two examples that the susceptibility can be expressed in terms of the scaling variables corre- 
sponding to the two fixed points governing the HT flows. We have given convincing evidence that the expansions of 
these variables have overlapping domains of convergence. Several interesting questions remain to be discussed. 

We have not discussed in detail the initial approach of the unstable fixed point. We have just shown that it can in 
principle be incorporated by calculating the subleading corrections in 7/2, 2/3 etc... This study depends on the details of 
the initial measure. A particularly interesting set of initial measures are the Landau-Ginzburg models in the vicinity 
of the Gaussian fixed point. We can use the usual perturbation theory in the quartic (or higher orders) coupling 
constant to construct the scaling variables associated with the Gaussian fixed point following the procedure described 
above, and interpolate between the Gaussian fixed point and the unstable fixed point. 

The question of small denominators and resonances have been circumvented by using only quantities which are free 
of this problem (namely yi and yi ) in our calculation of the susceptibility. The effects of the small denominators of 
the other scaling variables is a complex topic presently under study Q . 

We emphasize that the complete knowledge of the scaling functions and their inverse provides analytical expressions 
for all the thermodynamic quantities at any volume (see Eq. (4.16) in the simplified model). In practical calculations, 
one is naturally led to combine different expansions and it is thus important to use sets of interactions which are 
compatible, as done in Ref. B for gauge theories. 
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